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In contrast to the classical world, an unknown quantum state cannot be cloned ideally, as stated by the 
no-cloning theorem. However, it is expected that approximate or probabilistic quantum cloning will be 
necessary for different applications, and thus various quantum cloning machines have been designed. Phase 
quantum cloning is of particular interest because it can be used to attack the Bennett-Brassard 1984 (BB84) 
states used in quantum key distribution for secure communications. Here, we report the first 
room-temperature implementation of quantum phase cloning with a controllable phase in a solid-state 
system: the nitrogen- vacancy centre of a nanodiamond. The phase doner works well for all qubits located on 
the equator of the Bloch sphere. The phase is controlled and can be measured with high accuracy, and the 
experimental results are consistent with theoretical expectations. This experiment provides a basis for 
phase-controllable quantum information devices. 

The no-cloning theorem 1,2 is fundamental in quantum mechanics and quantum information processing. To 
clone quantum states, various quantum cloning machines, which can clone quantum states either approxi- 
mately 3 " or probabilistically 9 , have been proposed and implemented in different physical systems 1012 , see 
two reviews 1314 . In addition to exploring the capabilities of cloning processes in the quantum world, one direct 
application of quantum cloning is its role in analysing the security of quantum key distributions (QKDs), which 
are secure under all conditions because of the no-cloning theorem. The well-accepted Bennett-Brassard 1984 
(BB84) quantum key distribution protocol 15 uses four quantum states: to transmit information through a public 
quantum channel so that secure keys can be generated. A quantum strategy that can be adopted by an eaves- 
dropper is to quantum clone the transmitted BB84 states. Quantum cloning the BB84 states is equivalent to 
cloning the four states (|0) + e" p \ 1))/ s/l, where <p = jn/2,j = 0, 1, 2, 3, which we also refer to as BB84 states in this 
work. Note that all four of these states are located on the equator of the Bloch sphere. 

In the general quantum cloning process, a unitary transformation acts on a system composed of the BB84 input 
states, the blank state and the ancillary state, to obtain two optimal output copies. The general form can be written 
as U\00R) = a\00A) + b\01B) + c|10C) + d|UD), U]10R) = e|00£) + j\01F) + g\WG) + h\UH). Because the 
BB84 states are used randomly, an attack by quantum cloning should be symmetric in those four states. In general, 
we use the fidelity between the input state and the output copies to quantify the similarities between the input state 
and the output copies, which can be used as a measure of quality for a cloning machine. If two copies are 
equivalent, optimal cloning of the BB84 states is realised by the transformations 1 00) — » 1 00) and 
1 10) ->(|0l) + 1 W})/\/2. Here, the ancillary state is not necessarily needed; elimination of the ancillary state 

can conserve limited qubit resources. The corresponding optimal fidelity is F = — I — = Interestingly, we note 

2 v 8 

that both the cloning transformations and the fidelity presented above are optimal if we clone all qubits located on 
the equator of the Bloch sphere with the form (|0) +e" p \ 1)) l\2; here, the phase (p is arbitrary 71618 . Thus, the 
phase quantum cloning machine is optimal for cloning BB84 states; i.e., the phase quantum cloning machine is 
equivalent to the optimal cloning machine for BB84 states. 

For the previous experiment in a solid-state system 1 ', only two different phases were precisely prepared, and 
only a fixed phase parameter was measured. The arbitrary phase was realised through free time evolution. In this 
article, all of the different phases in our reported experiment are precisely prepared and controlled; this approach 
ensures that the phase cloning machine is successfully implemented, and the measurement is performed by 
standard state tomography, which is easier to implement. Our experiment is implemented in the nitrogen- 
vacancy (NV) defect centre of diamond at room temperature. This solid-state system is one of the most promising 
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candidates for quantum information processing (QIP) 20,21 . Many 
coherent control and manipulation processes have been performed 
with this system 22-34 . It is more challenging to implement the experi- 
ment in nanodiamond than in a bulk sample because of the coher- 
ence time. However, nanodiamond shows more potential for further 
integration and applications. Compared with a bulk sample, the NV 
centre in nanodiamond has many outstanding characteristics. 
Because the NV centre in a nanodiamond is sub-wavelength in size, 
reflection effects at the surface are avoided; thus, the fluorescence 
collection efficiency is much higher in nanodiamonds than in bulk 
materials 35 . In addition, with the aid of microfabrication techniques, 
the NV centres in nanodiamond can be manipulated and coupled to 
other systems, such as nanowires 36,37 , photonic crystals 38 , wave- 
guides 39,40 , microspheres 41 and optical fibres 42 . As a result, nanodia- 
monds play an important role in many applications, such as 
nanoscale imaging magnetometry, in which a nanodiamond is 
attached to an atomic force microscope (AFM) tip to measure the 
magnetic field 27,43 . The use of nanodiamonds has been demonstrated 
in plasmon-enhanced single photon emission through the assembly 
of individual nanodiamonds and gold nanoparticles 44 . Fluorescent 
nanodiamonds also offer many advantages in biophysics because of 
their brightness, photostability and biocompatibility 45 . 

Results 

Experimental system. Our experiment is performed in a 
nanodiamond sample, and an AFM scanning image of the sample 
is shown in Figure 1(a). The bright spots are nanocrystals with an 
average diameter of approximately 50 nm. The laser scanning image 
of the NV centre investigated is shown in Figure 1(b). The typical 
fluorescence intensity of nanodiamond is 10 s counts/s, which is 
approximately twice that of the bulk sample. Nanodiamond has 
many advantages but also poses new challenges for coherent 
control. The NV centre in nanodiamond typically has a short 
dephasing time, which is caused by the complex spin environment 



(rather than the 13 C spin bath in a bulk sample) 46 . The typical of 
the NV centre in nanodiamond is only hundreds of ns; Figure 1(c) 
shows the Ramsey signal of the investigated NV with ~ 700 ns. 
Thus, high-speed manipulation of the electron spin state is crucial for 
applications of the spin state to QIP, especially in complex systems 
with large microwave losses. In our experiment, a carefully calibrated 
coplanar waveguide (CPW) antenna is employed to deliver the 
microwave (MW) pulse to the sample [Fig. 1(b) inset]. The 
nanodiamonds are spin-coated on a CPW with a quartz substrate. 
The driven Rabi oscillation frequency reaches more than 50 MHz, 
which corresponds to a period of 20 ns. This duration is far shorter 
than the dephasing time of most NV centres in nanodiamonds. Thus, 
more state manipulations can be implemented before decoherence, 
Figure 1(c). 

The energy level structure of the electron spin of the NV is shown 
in Figure 1(d) with the corresponding optically detected magnetic 
resonance (ODMR) spectra. Levels m s = + 1 and m s = — 1 are split 
by approximately 250 MHz via an external magnetic field. The 
quantum states corresponding to m s = 0, m s = + 1 and m s = — 1 
are denoted by |0), | + 1) and | — 1). Laser and microwave pulses are 
used to manipulate the electron spin of the NV. Two types of micro- 
wave frequencies are used in the experiment: MW1 and MW2. The 
MW1 frequency is resonant with the transition |0)<->| + l), and 
MW2 is resonant with |0)<->| — 1). The MW1 signal is split into 
two parts: MW1-P and MW1-M, for state preparation and measure- 
ment, respectively. See the Methods section for more information 
about the experimental system. 

Experimental implementation. The implementation of this phase 
cloning machine includes preparation of the BB84 states, the cloning 
process and readout of the output. Three energy levels of the electron 
spin of the NV are used to implement this scheme. The logical BB84 
states are encoded by physical spin states as follows: j + 1) ~ |00) ( , 
|0) ~ 1 10), and | - 1) ~ |0l) ( 19 . The electron spin state of the NV can 
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Figure 1 | Nanodiamond samples and the energy level structure of the NV centre, a. AFM scanning image (30 jitn * 30 ;im) of the nanodiamond 
samples on CPW. The bright spots are nanocrystals, and the two straight bands are the metal edges of the CPW. b. The NV centre investigated is shown as a 
bright spot marked with a circle in the two-dimension laser scanning image, in which the CPW margins are also visible. The inset shows the CPW antenna 
fabricated on a quartz substrate, c. Typical Ramsey signal showing that TJ of the NV is approximately 700 ns. With high-frequency Rabi oscillations 
(inset), more manipulations can be implemented before decoherence. d. The three energy levels involved in the experiment and their corresponding 
ODMR spectra. 
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be polarised through laser excitation and can be read from the 
fluorescence intensity. A microwave pulse with duration 0 = 2nt/T 
(where t is the MW period and T is the Rabi period) can rotate the 
electron spin by an angle 0, around an axis in the equatorial plane of 
the Bloch sphere. The angle between the axis of rotation and the x- 
axis is determined by the phase of the microwave pulse. Thus, a 

0-0 

superposition state ji//) = cos- 10) +e" / ' sin- 1 1) can be prepared 

with a MW pulse of duration 0 and phase <p, Figure 2(c). 

As a result, the experiment can be performed by following the 
control sequence shown in Figure 2(a). First, the electron spin is 
polarised into state |0) with a 3- (is laser pulse. A MW1-P pulse with 
a duration of n/2 prepares the superposition state (|0) +e" p \ + 1)) 
/ \fl, where the variable <p is determined by the phase of MW1-P. 
This step combines the preparation of the input state for the first 
qubit and the initialisation of the blank state of the second qubit for 
the logical qubits. We note that the prepared input state actually 

corresponds to the form (|l) +e" p \0})/ \fl for the logical qubit; this 

form is equivalent to the normal form (|0) +e" p \ 1)) I \fl when (p 

ranges over all four values in the experiment. In the cloning process, 

an input state (|0) +e" p \ + l))/\/2 is transferred into an output state 

-|0)+ -| -1) + —=e i ' p \ J rl) byaMW2 pulse with a duration of nl 
2 2 y 2 

2. To measure the state phase and probability, a MW1-M or MW2 
pulse of variable duration is applied to induce Rabi oscillations 



between |0)<->| + l) or |0)<->| — 1), see the Methods section and 
Figure 2(d) for details. Finally, the spin state is read by photon count- 
ing under pulsed laser excitation. The readout signal is followed by a 
reference. The entire process is repeated to obtain the statistics of the 
measurements. 

Phase control and measurement. The accuracy of the control of the 
microwave phase and the measurement of the state phase is verified 
before cloning. Two mechanical phase shifters, with an adjustment 
range of 180°/GHz, are used to control the microwave phase. In 
addition to the readings from the phase shifters, the relative phase 
difference between the two microwaves can be estimated with an in- 
house designed device to measure microwave phase differences. A 
motor-driven phase shifter controlled by a computer modulates the 
MW1-M signal. During the experiment, the shifter automatically 
switches between two specific positions (X and Y). There is a 90° 
phase difference between microwaves at these positions of the shifter. 
The corresponding microwave phases are denoted Phase X (0°) and 
Phase Y (90°), which provide references for all of the following state 
phase measurements. 

The other phase shifter is manually controlled to modulate the 
MW1-P signal. Microwaves with various phases can be obtained by 
adjusting the manual phase shifter. Because a quantum state with a 
specific phase can be prepared using a microwave pulse with a con- 
sistent phase, three states with different phases are tested first. These 
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Figure 2 | Scheme for quantum phase cloning, a. Control sequence for the laser, microwave pulses and photon counter. After polarisation with the laser 
pulse, a n/2 MW1-P pulse prepares the input state. The quantum cloning process is performed by a n/2 MW2 pulse. MW1-M or MW2 pulses of various 
durations are used to implement the Rabi oscillations for the measurement, b. Phase measurement for three test states. The blue columns show the 
experimental results, and the red columns display the expectation values. The error bars are the standard deviations from seven measurements, c. Bloch 

6 6 

sphere representation of the superposition state \ijr) =cos- |0) + e" ? sin - |l) and the Bloch vector P and itsX, Y and Z components P„ P y and P z . d. Rabi 
oscillations are used for the measurement. The blue lines are the standard Rabi oscillations, which are used for scaling. The scatter points show the 
experimental data, and the green lines represent the fitting results. The upper panel shows the probability measurements. The orange arrow indicates the 
probability of 10), and the green arrow indicates the probability of I — 1) (Rabi C) or I + 1) (Rabi S). The lower panel shows the phase measurements. The 
dark yellow line denotes the middle of the Rabi X or Rabi Y oscillations. The orange arrow indicates the relative distance corresponding to the observation 
ofP* orP*. 
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state phases are measured over two full Rabi oscillations by MW1-M 
with X and Y phases; see the Methods section and the lower panel of 
Figure 2(d) for details. The experimental values are 40.4° ± 0.8°, 
49.2° ± 1.2° and 60.5° ± 0.8° (standard deviation of seven measure- 
ments), whereas the expected values are 40°, 50° and 60°, Figure 2(b). 
These results demonstrate the high accuracy of the control and the 
measurement; the small standard deviations also show good phase 
stability for repetition at long times. 

Experimental results. The cloning experiment is implemented with 
four different input states that are distributed uniformly on the 
equator of the Bloch sphere, see Figure 3. Coincidentally, these 
states can serve as BB84 states in QKD. 

First, an input state <t>i of the form ( j 0) + e" p | 1 ) ) / \fl is prepared 

with the polarisation laser and a 7C/2 MW1-P pulse. The state phase is 
measured five times, and the experimental data points are shown as 
empty green squares in Figure 3. The average phase before cloning, 
44.6°, is shown by the blue line with solid triangles. The standard 
deviation of the five measurements is ±2.4°. In the subsequent clon- 
ing stage, the input state is converted into the -|0)-| — | — 1) 

1 2 2 

H — = e" p \ + 1) state with a nl2 MW2 pulse. The phase measurement 
\/2 

results for the output state are shown as empty red circles. The 
average phase after cloning, 50.4°, is represented by the magenta line 
with solid triangles, and the standard deviation is ± 1.9°. The phase 
change for the cloning process is 5.8°; this change is a reasonable 
value based on the magnitude of the deviation in the measurement. 
Thus, the phase information is well preserved during the cloning 
process. 

To confirm the cloning fidelity, we are concerned with not only the 
state phase but also the probabilities of |0), | + 1) and | — 1). The 
probability measurement is also performed five times; the detailed 
process is described in the Methods section and shown in the upper 
panel of Figure 2(d). The probabilities of |0) and | + 1) can be esti- 
mated in a similar manner from the Rabi oscillation between the two 
states. The measured results are (25.4 ± 1.3)% for |0) and (51.1 ± 
2.1)% for | + 1). The probabilities of |0) and j — 1) can be estimated by 
the Rabi oscillation between the two states. The measured results are 
(25.9 ± 0.8)% for |0) and (24.7 ± 1.5)% for |-1). Summarised 
together, the relative probabilities of |0), | + 1) and | — 1) can be deter- 
mined and are shown in inset "a" of Figure 3. The probability values 
are 25.22% for |0) (yellow sector), 50.44% for | + 1) (red sector) and 
50.44% for | + 1) (green sector). For the output state with the form 

1 1°) + \ I - !) + -r eiV \ + !>> the expected values are 25%, 50% and 

2 2 y2 

25%. Based on the magnitude of the deviation of the measure- 
ment, the experimental results are consistent with the expected 
values. 

These experimental procedures are repeated with three other 
states: <5 2 , (b^and®^ States with different phases are prepared 
through adjustments of the manual phase shifter, and the state 
phases are measured before cloning. The measured results are 
128.2° ± 1.6°, -47.7° ± 1.7° and -138.4° ± 2°. The state phase 
measurement is performed again after the cloning process. The phase 
values are 124.3° ± 1.7°, -48.3° ± 3.3° and - 140.3° ± 2.9°, respect- 
ively. The errors are the standard deviation of five measurements, 
and these results are shown in Figure 3. The phase change values 
during the cloning process are —3.9°, —0.6° and —1.9°. The phase 
information is well preserved during the cloning process. The prob- 
abilities of states <P 2 , Q>jandQ> 4 are also confirmed by measurements. 
The relative probabilities of |0), | + l) and |-1) are 23.26%, 52.49% 
and 24.24% for state 0 2 (Fig. 4, inset b), 25.7%, 50.05% and 24.25% 
for state <D 3 (Fig. 4, inset c) and 24.95%, 49.23% and 25.82% for state 
<1>4 (Fig. 4, inset d). All of the results are highly consistent with the 
expected values. 
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Figure 3 | Measured results for the quantum state phase and relative 
probability. The scatter points show the experimental data, and the lines 
with the solid triangles display the averaged values. The blue lines represent 
the measured phases CD;, <J> 2 , ft>^and<S> 4 for four different input states, and 
the magenta lines show the measured phases after the cloning process. 
Insets a, b, c and d show the relative probability results after cloning for 
states 3>i, ® 2 > ( S>iand<S> 4 , respectively. The yellow sectors show the 
probability of 1 0), red shows I + 1), and green shows 1 — 1). 

Density matrices and fidelities. We can denote the physical output 
state in the forma|0) +/ie i( °| + 1) +y| - 1) and |a| 2 + \p\ 2 + |y| 2 = 1, 
because the relative phase of | — 1) is fixed at 0° in all of the 
experiments. From the results of the state phase and probability 
measurements, the physical output states and their density 
matrices can be determined. The matrix elements of states Oj, <t> 2 > 
Q) 3 andQ> 4 are shown in Figure 4. The diagonal elements have a one- 
to-one correspondence between the matrix and the columns in the 
figure. The yellow columns indicate the elements p 0 0 , red indicates 
P+i,+i) an d green indicates Because the density matrix is 

conjugated, the real and imaginary parts of the off-diagonal 
elements are shown together in Figure 4. The blue columns 
indicate the real part (p 0 ,+i + p+i j0 )/2, and the cyan columns 
indicate the imaginary part (p + i )0 — p 0 ,+i)/2i. Matrix elements 
that are not involved in the experiment are shown as grey columns 
with calculated heights. 

Decoding the logic output states from the physical states, we can 
calculate the fidelities between the input states and their copies after 
cloning. By tracing out one qubit degree of freedom from the density 
matrix of the output state after decoding, we obtain the two copies 
(p x and p 2 ) of the input state Mi/r) = (|0) +e ilf \ 1)) j . In mea- 
suring the phase of the input state before cloning, we note that the 
sum of P x and P y 2 is less than 1, (P z = 0, see Fig. 2(c) for information 
about P x , P y , P z ). Phase amplitude (P x 2 + P y 2 ) damping occurs during 
preparation and measurement 47 . We scale P x 2 + P y 2 to 1 and com- 
pensate the results for P* and P* after cloning based on the damping 
rate. For the compensated density matrix p r , the fidelities of the two 
copies are F\ = {\l/\p\ Ji/^andFi; = (\j>\p\ \\ lA)- The calculated results 
are F\ = 78.5%andF L 2 = 78.0% for state <S> U F[ = 79.7%and 
F L 2 = 80.3% for state <P 2 , F[ = 84.6%andF c 2 = 83.0% for state O s and 
F[ = 79. 0%andF 2 = 79.5% for state < 5 4 . The compensated experi- 
mental fidelities are 78.3% for state 80.0% for state 3> 2> 83.8% 
for state {D 3 and 79.2% for state <J> 4 . The average experimental fidelity 
is 80.3%. 
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Figure 4 | Results for density matrices and phase cloning fidelity, a, b, c and d illustrate the results for states ®i, <J> 2 , <J> 3 a«d<I>4, respectively. 
The yellow columns correspond to the diagonal elements of the density matrix and indicate p 00 , red indicates and green indicates P-i-i- The real 

and imaginary parts of the off-diagonal elements are shown together. The blue columns indicate the real part (po,+ i + p+i,o)/2> and the cyan columns 
indicate the imaginary part (p+i_o — Po,+i)/2z- The matrix elements that are not involved are illustrated by grey columns with calculated heights, 
e. Fidelity results for states <S> 1 , <J> 2 , <S>3and$> 4 and their average values. The blue and green columns show the fidelities of the two output copies; 
the averages for these copies are shown with magenta columns. The red dash-dot line represents the expected optimal cloning fidelity, 85.4%, 
and the orange line displays the universal cloning bound, 83.3%. 



The extracted information for the relative probabilities and phase 
(p is sufficient for phase cloning, so we can calculate the phase fidelity 
separately. This calculation includes the additional assumption that 
phase decoherence is negligible in the experiment. If phase damping 
is ignored, the fidelities of the two output copies are F l = (i/'|Pi|i/') 
and F 2 = {iplPiW- The calculated results are F l = 84.9%andF 2 = 
84.3% for state <Di, F 1 = 84£%andF 2 = 85.3% for state 4> 2 , F i = 
85.9%andF 2 = 84.8% for state <TJ 3 and F 1 = 85.0%andF 2 = 85.6% for 
state €> 4 . The phase fidelity results are 84.6% for state Oj, 85.0% for 
state <1> 2 , 85.4% for state €> 3 and 85.3% for state tf> 4 . The average phase 
fidelity, 85.1%, is very close to the theoretical fidelity, 85.4%. These 
results are shown in Figure 4(e). Thus, phase cloning is successfully 
implemented in the experiment. 

Discussion 

In this work, we have demonstrated the first room-temperature 
implementation of a quantum phase cloning machine with control- 
lable phase in a solid-state system composed of nanodiamond. 
Nanodiamond has great advantages for future applications; however, 
its short dephasing time brings new challenges. A carefully calibrated 
CPW antenna was employed to improve the microwave delivery 
efficiency. For high-frequency Rabi oscillations, more manipulations 
can be implemented before decoherence. The microwave phase 
is precisely controlled in our experiment, and the phase of the 
quantum state can also be measured with high accuracy. With these 



technical developments and preparations, we implemented optimal 
phase-covariant quantum cloning in our NV nanodiamond system. 
The experimental results were measured by standard state tomo- 
graphy methods, which are more precisely. The phase information 
is well preserved during the cloning process. Decoherence from 
phase damping causes the average experimental fidelity of cloning 
to be 80.3%, which is 94% of the theoretical value. The fidelity can be 
improved by a NV with longer dephasing time, faster manipulation 
and dynamic decoupling techniques 48 51 . For a phase cloning ma- 
chine, under the assumption that phase damping in the experiment is 
negligible, we can use only the relative probabilities and the phase <p 
to analyse the cloning efficiency of the information. In this way, the 
calculated phase fidelity is 85.1% on average; this value is close to the 
theoretical value, 85.4%, and exceeds the universal quantum cloning 
bound, 83.3%, Figure 4(e). The object states cloned in the experiment 
can serve as BB84 states in QKDs. Thus, this high-fidelity cloning 
may play an important role in analysing the security of QKDs. The 
drawback of this quantum cloning experiment is that the two qubits 
are encoded in a one-spin system. One possible method for solving 
this problem is to use both the electron spin and the nuclear spin of 
the NV centre to encode the two qubits; this encoding is still chal- 
lenging because the two spins have different coherence times and an 
accurate readout is difficult. The full control and precise mea- 
surement of the quantum phase can be used as a foundation for 
further applications in QIP. The coherent manipulation and 
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implementation of the cloning machine in nanodiamond also dem- 
onstrate prospects for future scalable quantum information devices. 

Methods 

Experimental setup. A single NV centre can be found and located with our home- 
built laser confocal scanning microscope system. The fluorescence emitted from the 
NV centre is collected with the optical system and is then detected by a single photon 
counting module connected to the counter terminal of a multifunction data 
acquisition device (NI USB 621 1). To obtain the laser pulses for state polarisation and 
for readout, an acoustic optical modulator (AOM) modulates a 532-nm continuous- 
wave laser. The microwave signals are switched by mixer (Minicircuit) pairs for a 
quick response (nearly real-time) and for a better on-off contrast ratio. The signals 
that control the mixers are provided by pattern/pulse generators (Agilent 81110A), 
which run in triggered pulse mode with high time resolution (better than 0.1 ns). 

State phase measurement. The Bloch vector for one quantum state can be rotated 90° 
about the x- ory-axis by a nil microwave pulse with phase X or Y 52 ' 53 , [Fig. 2(c)]. The 
X or Y component (P x , P y ) of the Bloch vector P is converted to the Z component (P z ), 
which can be measured through photon counting. Full Rabi oscillation using a 
microwave with phase X or Y, Rabi X or Rabi Y, is implemented in the experiment. 
The red circles filled with yellow in Figure 2(d) show the experimental data. Rabi X 
and Rabi Y are scaled by the standard Rabi oscillation from the pure 1 0) state (shown 
by the blue line). In the lower panel of Figure 2(d), the orange arrow indicates the 
relative distance between the point at the n!2 position of the fitting line for the Rabi X 
or Rabi Y data (green) and the line at the middle of the oscillation (dark yellow). Based 
on the sign (positive or negative), the distance indicates the observation of P x or P y in 
the Bloch sphere representation. After the cloning process, all three energy levels are 
involved, and the value is denoted as P* or P* . After the measurement of Rabi X and Y, 
the state phase can be obtained from the value of arctan P* j ' P x . 

State probability measurements. The fluorescence intensity under pulsed laser 
excitation indicates the probability of 1 0), which is easily defined by the initial point of 
the Rabi oscillation. The probability measurement for | + 1) or | — 1) is indirect, 
however, and is performed by rotating the Bloch vector and reversing the Z 
component with a n pulse from MW1 or MW2 [Fig. 2(c)] . Two full Rabi oscillations, 
Rabi C and Rabi S, are implemented in our experiment. Rabi C with MW2 is used to 
estimate the probability of |0) and | — 1); Rabi S with MW1-M is used to estimate |0) 
and | + 1). Based on prior state phase measurement results, the microwave phase of 
Rabi S can be adjusted to be similar to that of the prepared state. In the actual 
experiment, when the phase difference is too large, the fidelity of the measurement 
process may decrease. In the worse-case scenario, for example, if the two phases 
are nearly opposite, the fidelity may be less than 80%. As a result, phase 
adjustment based on prior phase measurement results is beneficial for high 
measurement fidelity. If scaled by the standard Rabi oscillation of the pure 
|0) state, the initial height of Rabi C or Rabi S gives the probability of |0), 
[Fig. 2(d), orange arrow in the upper panel] . The height of the fitting line at n gives the 
probability of | — 1) or | + 1), [Fig. 2(d), green arrow in the upper panel]. 
For the probability of the |0) state, the demonstrated values are the average of the 
Rabi C and Rabi S measurements. 
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